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Abstract 

We consider a model for systems perturbed by dichotomous noise, in which the 
hazard rate function of a random lifetime is subject to additive time-alternating per- 
turbations described by the telegraph process. This leads us to define a real-valued 
continuous-time stochastic process of alternating type expressed in terms of the inte- 
grated telegraph process for which we obtain the probability distribution, mean and 
variance. An application to survival analysis and reliability data sets based on confi- 
dence bands for estimated hazard rate functions is also provided. 
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1 Introduction 

Large attention has been given recently in the physical and mathematical literature to 
stochastic systems perturbed by dichotomous noise, such as those related to Brownian mo- 
tion or Brownian motors (see, for instance, Bena et al. [5], [1] and Porra et al. [H]). The 
relevance of Markovian dichotomous noise in biological systems has been also pointed out. 
See, for instance, Laing and Longtin [14] where the benefic role of noise is investigated 
within the context of mathematical neuroscience. In this paper we aim to study the effect 
of dichotomous noise to the hazard rate function of random lifetimes. 

It is well known that in many fields related to survival analysis and reliability theory 
the study of hazard rates plays a very important role (see the classical book of Barlow and 
Proschan [3] , for instance) . Specifically, here we are interested to hazard rate functions that 
are realizations of stochastic processes. Such types of hazard rates arise in some models 
that were recently proposed to describe doubly random phenomena such as lifetimes of 
devices operating in random environments. See, for instance, Kebir [13) and the included 
references for a model in which the hazard rate is a functional of a stochastic process that 
describes the random variability of the environment. A similar model is treated in Di 
Crescenzo and Pellerey [9] , where the hazard function is the realization of a non-decreasing 
stochastic process with independent increments. Other models characterized by stochastic 
hazard rate functions were also discussed by Ezov and Jadrenko [llj and by Yadrenko and 
Zhegrh |28j . Furthermore, a stochastic model in which a constant hazard rate (describing 
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the fractional dissolution rate of a drug potion) is corrupted by a white noise was recently 
analyzed by Lansky and Lanska [TB] and by Lansky and Weiss [T7] within the theory of 
drug dissolution. Other applications of stochastic hazard rates were proposed in the actuarial 
literature and in studies of population biology and human aging (see Milevsky and Promislow 
[T5] and references therein). All these stochastic models are suitable to describe phenomena 
characterized by an intrinsic randomness, which is evidenced at the hazard rate level. 

Along the line traced by the above mentioned papers, herewith we propose a new model 
to describe systems characterized by failure rates which are subject to random perturbations 
expressed as dichotomous noise. We aim to discuss a model that incorporates a more realistic 
kind of noise characterized by a constant intensity, and that at the same time presents the 
feature of noise with alternating behavior that is frequently found in biological systems. 
For instance, see the papers by Buonocore et al. [7j and Di Crescenzo et al. [8], dealing 
with neuronal dynamics subject to alternating inputs. Stochastic models characterized by 
alternating behaviour are also used to describe the cell motility on one-dimensional and two- 
dimensional state-spaces (see, for instance, Lutscher [TS], Stadje [55], and references 
therein) . 

Suppose that (non-negative) lifetime T has distribution function 
(1) F(t) = 1 -expj-J^ r(s)ds| , t > 0, 

so that its hazard rate function is 

/ s P(T < t + At I T > t) 
(2 r(t)= lim 1 = \ 1 '-, t > 0. 

w W At^0+ At 

We assume that r(t) > c, where c is a positive constant, and that r(t) is subject to ad- 
ditive time-alternating perturbations described by the well-known telegraph process (see, 
for instance, Orsingher \2Q\ or Beghin et al. [1]). This assumption consists of applying the 
substitution 

r(t) r{t) + V{0){-l) N(t \ t>0 

in the right-hand-side of (P), where V(0) S {— c, c} and {N(t);t > 0} is a Poisson process 
independent of T. We are thus led to defining a real-valued continuous-time stochastic 
process {X(t); t > 0}, where X(t) is expressed in terms of the integrated telegraph process: 

(3) A(t) = l-exp|-^*[r( S ) + ^(0)(-l) Ar ( s )]d S j, t > 0. 

Here, we present some basic results on process (j3|). First of all, in Section [2] we obtain 
the moment generating function of the integrated telegraph process, which is involved in the 
expression of mean and variance of X(t). Such moments and the probability distribution 
of X(t) are given in Section [3] This distribution has a discrete as well as an absolutely 
continuous component. Finally, in Section [4] we use an asymptotic confidence band for 
estimated hazard rates in two case-studies where our model provides adequate fit to two 
data sets taken from the survival analysis and reliability theory literature. 

We point out that the results obtained in this paper can be extended to the case where 
there are two different values V^O) € {— C2,ci} and N(t) is an alternating Poisson process 
characterized by two rates Ai and A2. Finally, we point out that stochastic processes ob- 
tained by transformations of the integrated telegraph process are of interest in stochastic 
modelling literature. For example, Di Crescenzo and Pellerey [10] use a geometric telegraph 
process to model the price dynamics of risky assets. 



2 Moment generating function of integrated telegraph 
process 

Let {W(t);t > 0} denote the well-known integrated telegraph process, defined by 



W(t)=V(Q) / (-!)*«<!», <>0, 



(4) 



where {N(t); t > 0} is a homogeneous Poisson process with intensity A and V(0) is a random 
variable independent from N(t) such that 



(5) 



P{V(0) = ±c} 



2' 



with c > 0. The distribution of (0| is characterized by discrete components concentrated at 
±ct and by an absolutely component on (—ct,ct). Orsingher [2D], shows that the discrete 
component satisfies 

(6) P{W{t)=ct}=P{W(t) = -ct} = -e- Xt , t>0, 

and that, for t > 0, ~ct < x < ct, the continuous component has density 

A Jo 0v / c 2 t 2 -.x 2 ) + ^J (^-Vc^~:- 



(7) °- Fm t)<*} = ±e-» 



where 
(8) 



Jo(*)=£ 



(x/2) 



2 k- 



denotes the modified Bessel function of order 0. 

Proposition 2.1 For all s G M and i > £/ie moment generating function of W{t) is 
(9) M(s,t):=E 



c 

-At 



W(t) 

cosh ^a/a 2 



A 



\/A 2 + s 2 c 2 



sinh 



+ s 2 c 2 



Proof. From Eqs. © and Q it follows 

„-Ai 



(10) M(s,t) 



(c sct + e- sct ) 



2c 
-xt r 



A 



AJ ( -VcW 



2c 



XQ(s,t) + —Q(s,t) 



dt 

s € R, t > 



d.z: 



where we have set 



QM) := 



y e sx I (^VcH 2 - xA dx. 



Making use of Eq. (25), from Orsingher [20] we obtain: 



X 2 Q(s,t). 



A two-fold integration by parts shows that 



d 



^Q{ s ,t) = (\ 2 + S 2 c 2 )Q{s,t). 



Solving this equation with initial conditions Q(s,0) = and — Q(s, t) | 4 _q = 2c we have: 



Q(a,t) 



VA 2 + s 2 c 2 



S S R, t > 0. 



Using this formula in the right-hand-side of (|10p. expression ([9]) finally follows. 



□ 



It should be noticed that © could also be obtained from the initial-value problem for 
the telegraph equation 

( a 2 n B 2 d 2 

—rp + 2X—p = C ■ TT-yP 



< 



dt 21 ' 1 ~"dt* 

p(x,0) = 5(x) 



dx 2 



I dt p{x>t) 



= 0, 



where S(x) is the Dirac delta function. Indeed, M is solution of 

f & M + 2A = s 2 c 2 M 
at z dt 



M(s,0) = 1 



= 0. 



3 Results 

Let T be a non-negative random variable, with absolutely continuous distribution func- 
tion {J), hazard rate function ([2]), and support 

S:=[0,£), where £ = sup{t > : F(t) < 1} G (0, +oo]. 

Throughout the paper we will assume that r(t) > c for all t £ S, where c is a positive 
constant. Hence, 

F(0) = 0, F(t) > 1 - c~ ct for all t > 0, 

so that 

T< st Exp(c), 

where '< s t' denotes the usual stochastic order (see, for instance, Shaked and Shantiku- 
mar |23j ) and Exp(c) denotes an exponentially distributed random variable with mean c _1 . 
Let us now consider the stochastic process {X(t); t £ S}, with X(0) = and 



(11) 



X(t) = 1-F(t)e- W ^, teS, 
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Figure 1: (a) Two sample-paths of process @ for c = 2 and A = 15; (b) the corresponding sample- 
paths of fTTj), with F(t) and r(t) defined in (fT2]) and ([13]) . respectively, with a = 15 and /3 = 0.001; 
the solid line shows 



where W(t) has been defined in (j4|) and where F(i) = 1 — F(t) denotes the survival function 
of T. Note that process (fTTj) is equivalent to process Hence, for any fixed t, X(t) 
describes the probability that a random lifetime is not larger than t. The system hazard 
rate is random and is given by r(t) + V(0)(-l) N ^\ where the distribution of 7(0) is specified 
in (0 . As pointed out in Section 1 , this hazard rate accounts for random perturbations on 
r(t) that occur according to a telegraph process. The sample-paths of X(t) are absolutely 
continuous distribution functions that approach |T]) as c goes to 0. As an example, Figure 
1 shows (a) two sample-paths of W(t), and (b) the corresponding sample-paths of X(t), 
characterized by distribution function 



(12) F(t) = 1 - exp 



-t i -—t 3 + -t 2 + (c + P)t 



t > 0, 



where 

(13) r(t) = at(t - l) 2 + c + /?, t>0, (a,/3>0). 

Next, we derive the probability distribution of X(t). Since p{— ct < W(t) < ct} = 1, 
t > 0, from Eq. (fTTj) we have: 

P{a(t) < X(t) < b(t)} = 1, t e 5, 

and 

a(t) < F(t) < b(t) for all t G (0, £), 

where 

(14) a(t) = 1 -F(i)e ct , = 1 - F{t) e~ c \ teS. 

From (|14p . we note that a(f) and b(t) are distribution functions with hazard functions r(t) — c 
and r(t) + c, respectively. Furthermore, there results: 

Hma(t) = 1 - e~", ]imb(t) = 1, 




Figure 2: Plots of D(t) in 3 different cases: (a) r(t) given by Eq. (fl3"Tl . with c = 1, A = a = 15 and 
/3 = 0.001; (b) r(t) = l + e*, t > 0, with c = A = 1; (c) r(i) = 1 + 3(1 - e _t )/(e i + e _t ), t > 0, with 
c = A = 1. 



where 

[r(s) — c] ds. 

We note that v = +oo when I < +oo, whereas v < +oo if £ = +oo. Moreover, by setting 
D(t) := - a(£), from ([Til) we obtain: 

D{t) = 2F(t) sinh(ct), t e <S, 

with -D(O) = and \im D(t) — e _l/ , with D(t) non-decreasing in t if 

r(t) < c coth(ci). 

For instance, D(t) is shown in Figure 2 for three different choices of r(t). The first choice 
shows a case in which D(t) is bimodal. In the third choice we have lim D(t) — 0.268, 

t — >-f oo 

whereas such limit vanishes in the first two choices. 

We shall now determine the probability distribution of X(t). Like as for the integrated 
telegraph process, the distribution of X(t) consists of a discrete component on a{t) and b(t) 
and an absolutely continuous component inside of (a(f),6(t)). 

Proposition 3.1 For all t G S the discrete component of the distribution of X(t) is: 
(15) p{X(t) = a(t)} = P {X(t) - b(t)} = 1 c 



-At 



Moreover, for all x € (a(t),6(t)J and t G S the continuous component has density 
(16) /M) = ^7^e 



-A / 



XIq f^^u(x,t)j + (^y/u(x,t) 



where 
(17) 



u{x, i) := c 2 t 2 - In 
= In 



F(t) 
1-x 
1 - a(t) 



1-x 



In 



1-x 



with a(t) and b(t) defined in (I14p . 

Proof. Due to (fTTj) . the distribution function of X(t) can be expressed as 

F(t) 



P{X(t) < x} = P <^ W(t) < In 



1-x 



t G 5, a(t) < x < 6(f). 



Hence, recalling © and 0, Eqs. {TSJl and JTSJ) follow. 



□ 



The expression of density p^|) is similar to that of the underlying integrated telegraph 
process. Unlike the latter process, however, the support of f(x,t) is bounded when t in- 
creases. Moreover, the expression for f(x, t) reflects the form of interval (a(t), b(t)). Indeed, 
from (117|) it is not hard to note that u(x,t) — when x — a(t) and x = bit), whereas 
u(x,t) > when x £ (a(t),b(t)). Furthermore, recalling (fl4|) and ([8|), for t € S we have: 



,. *, s Ae-( A+C )* / At 

lim f(x,t) = = H 

xia(*) 2cF(t) V 2 



lim f(x,t) 



Ae -(A-c)t 

2cF(t) 



1 



Xt 



Figure 3 shows for instance f(x,i) for three different choices of t. 
Proposition 3.2 For all t > mean and variance of X(t) are: 

(18) E[X(t)] = l-F(t)M(-l,t), 

(19) Var [X(t)\ = F 2 (t) {M(-2, t) - [M (-1, i)] 2 } 

where M(s,t) is given in (0). 
Proof. Recalling (|11|). we have 

E [*(«)] - 1 - F(t) E ^e _w/(t) 

from which Eq. (fT5|) immediately follows. Moreover, due to 

E [X 2 (t)] = E (l - 2F(t) e~ w ^ + F 2 (t) e~ 2W ^ 



= 1 



2F(i)E(c-^ t >) +F 2 (t)E(e~ 2W ^ 



we obtain: 



Var [X(t)} = 1 - 2F(t) E (e~ w ^ + F 2 (t) E (c 

-l + 2F(i) E (e- w W)- [e (e- ff «) 
! w{E(e- 2H/ (*))-[E (e-^>) 



2 



F*(t) 



F 



and thus Eq. J19] 



□ 
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Figure 3: Plot of density (|16|l for three choices of t and for x ranging within (0, 1); r(t) and all 
parameters are chosen as in case (a) of Figure 2, i.e. c — 1, X = a = 15, (3 = 0.001 and r(t) is given 
by Eq. ((15). 



Proposition 3.3 The mean ofX(t) is increasing int £ S, and is decreasing inc. Moreover, 



(20) 
(21) 



Proof. Using Eqs 

d_ 

dt 



lim E[X{t)] 
t^o 

lim E[X(t)] 
© and ©, 



0. 

1, 



lim Var[X(t)] = 0; 
lim Vax[X(t)] = 0. 



r F(t)M (-!,<)] 



r(i) + A - VA 2 + c 2 



A 



r(t) + X + VA 2 + c 2 



VA 2 + c 2 
A 



exp 



r(u)du - At + V A 2 + c 2 t 



VA 2 + c 2 



exp 



r(u)du - Xt - VA 2 + c 2 t 



Hence, recalling r(t) > c for all t € S we have i[F(i)Af(-l,i)] < 0. By virtue of (18 



the mean of X(t) is thus increasing in t. Moreover, from Eq. © it is not hard to see that 
M(s,t) is increasing in c > 0. Hence, E[-X"(t)] is decreasing in c due to ((T8| . Limits ((20)) 
easily follow from (fT8|) and ((T9|) . Recalling Eqs. (fT8| and Q we have: 



lim exp 



r(u)du - Xt + \/X 2 +c 2 t \ , 



where 



lim exp 

t->£ 



j-jf r(u)du - At + VA 2 + c 2 t 



(22) 



lim exp <j — I [r(u) — c] duf exp \ — j c + A — \/ X 2 + c 2 



du > = 0. 



The first of (j2Tj) thus follows. Moreover, due to Eqs. (fT9|) and (O, we have: 



limVar[X(t)] = -(l 



A 



x lim exp 

t->t 



VA 2 + 4c 2 



r(u) — c] du f exp 

2 



2c + A - ^A 2 + 4c 2 



du 



x lim exp /— 2 J [r(u) — c] du| exp |— 2 ^ c + A — \/ A 2 + 



du 



Finally, noting that 



lim exp 



and recalling Eq. 



f [r(u) - c] duj exp |- [2c + A - ^A 2 + 4c 2 ] dttj 



the second of (I2T1) holds 



□ 



Figure 4 shows mean and variance of X(t) for an example involving distribution func- 
tion (fl2|) and hazard rate (fT3|) . 

An immediate consequence of Proposition 13.31 is the following 

Proposition 3.4 Process X(t) converges in probability to I as t — > I. 

Finally, denoting by V(t) — V(0) (—l) N ^ the telegraph process that describes the per- 
turbating noise, making use of a well-known result (see, for instance, Theorem 3.4.4 of Ross 
[22] ) the following asymptotic probabilities follow: 



lim P {V(t) =c}= lim P{V(t) = -c} = - 

t — *+oo t — »+oc Z 



4 An application 

In this section we indicate a statistical procedure based on an asymptotic confidence band 
for estimated hazard rates, which is useful to assess the validity of the proposed stochastic 
model. 

Let Ti, T2, . . . , T n be iid absolutely continuous random variables describing a random 
sample of failure times having density f(t) and distribution function F(t). The usual kernel 
estimator of f(t) is 





Figure 4: (a) Mean and (b) variance of X(t); r(t) and all parameters are as in case (a) of Figure 
2, i.e. c = 1, A = a = 15, (1 = 0.001 and r(t) is given by Eq. JT3J. 



where k(-) is a bounded even density function, and h is the bandwidth. The corresponding 
estimator of F(t) is 

(24) -£>>*•- Sg*^ 

where K{t) = f ^ k{u) du. From (23]) and d24j) we can build up the following estimator for 
the hazard rate function: 

f(t) 



(25) 



f(t) 



l-f(i) 

Taking into account the asymptotic normality of l|23p and (1241) . with Var(/) = (n/i) _1 /C/, 
where /C = k 2 (t) dt, and Var(F) = 0(n _1 ), it follows that also r(t) is asymptotically 
normal, with variance Var(r) = (nh)^ 1 JCr 2 / f (see Hall et al. 12J). We can thus consider 
the following confidence bands for r(t), having nominal coverage 1 — 2a, 



(26) f-(t)=f(t) 



K 



nhf(t) 



1/2 



f(i)z Q , r+(t)=f(t) + 



K 



nhf(t) 



1/2 



where z a is the a-level critical point of the standard normal distribution. 

We note that a random hazard rate of type r(t) + V(0)(—1) N ^\ related to process ([3]), 
alternates randomly between functions r(t) — c and r(t) + c, thus being contained inside 
the strip r(t) ± c. This suggests the following procedure to test if the random sample 
(T 1 ,T 2 , . . . ,T n ) can be viewed as drawn from the (random) distribution function X{t): If, 
for a fixed a, the strip r(t) ± c is contained within a realization of the confidence band (f^Hj) 
obtained from observed data, the model defined by is defensible for the given data set. 




Figure 5: For the sample data of Application with r(t) = 0.0125, t > 0, a = 0.025, c = 0.0004, 
h — 6 and k(t) given in (|28p . it is shown: (a) the estimated density obtained by Eq. (|23p . (b) the 
estimated distribution function obtained by Eq. (|24[) . (c) the realization of the confidence band (|26|) . 
(d) the realization of the functions appearing in the left-hand-side (dashed line) and in the right- 
hand-side (solid line) of Eq. (gTJ). 

In other terms, with a (1 — 2a)-level confidence, we can adopt model X(t) if 

1/2 

r{t) z a - c for all t e (t {1 . n) , t( n _ 1:n) ), 

where denotes the value of the j-th order statistic in the observed random sample. We 
stress that, due to (j2"Tj) , this procedure is effective to reject "large" values of c. 

Application 4.1 — A case-study with constant hazard rate function. 

Let us consider the following set of 46 sample data given in Ahmad pQ : 

{13, 14, 19, 19, 20, 21, 23, 23, 25, 26, 26, 27, 27, 31, 32, 34, 34, 
37, 38, 38, 46, 46, 50, 53, 54, 57, 58, 59, 60, 65, 65, 66, 70, 85, 
90, 98, 102, 103, 110, 118, 124, 130, 136, 138, 141, 234}. 

These represent the survival times of certain patients in a melanoma study conducted by the 
Central Oncology Group (see Susarla and Van Ryzin [27] ). Figure [5^a) shows the estimated 
density obtained by use of (f2"3"| . where k{t) is the Epanechnikov kernel (see Silverman [2"4] ) 

(28) *(') = {iM 1_ £) 

l> otherwise. 

The corresponding estimated distribution function is given in Figure [5](b) . A bandwidth 
depending on the sample standard deviation was proposed in Azzalini [2], where it is also 



(27) 



\r(t)-r(t)\ < 



nhf(t) 



12 




Figure 6: For the sample data of Application 14.21 with r(t) given in J29), a = 0.025, c = 0.00025, 
h = 75 and k(t) given in (|28p . we show the same functions appearing in Figure [5] 



mentioned that these results are hardly affected if the standard deviation is estimated from 
the sample data. Hence, we have empirically used the fixed bandwidth h = 6; this choice has 
been also motivated by the need of obtaining a sufficiently smooth estimated density with 
small tails. We consider a baseline constant hazard rate function r(t) = 0.0125, t > 0. The 
realization of the confidence band (|2"rj|) obtained from the above mentioned data is plotted 
in Figure03c) for a — 0.025. With a 0.95-level confidence the strip r(t)±c, with c = 0.0004, 
falls within the confidence band, i.e. condition (|2~T|) is fulfilled (this is graphically shown 
in Fig. Od)). In conclusion, stochastic model §5§ is defensible for the observed data when 
r(t) = 0.0125, t > 0, and c < 0.0004. 

Application 4.2 A case-study with varying hazard rate function. 

Hereafter we analyze the following set of n = 86 sample data, taken from Table 1 of Langseth 
and Lindqvist [T5] : 

{220, 233, 234, 240, 265, 270, 273, 279, 285, 287, 294, 295, 300, 325, 328, 
333, 365, 368, 369, 381, 417, 418, 429, 460, 470, 474, 475, 476, 508, 522, 523, 
535, 542, 570, 580, 604, 612, 613, 614, 615, 634, 636, 637, 638, 651, 657, 660, 
666, 668, 680, 681, 684, 691, 693, 705, 717, 834, 837, 841, 843, 845, 875, 972, 
1037, 1084, 1091, 1109, 1117, 1197, 1258, 1269, 1297, 1309, 1322, 1346, 1349, 
1359, 1363, 1448, 1476, 1481, 1557, 1606, 1610, 1642, 1659}. 

They describes the service time of a single component in a reliability study. Figures [f^a) 
and[6^b) show the estimated density obtained from (|23[) . with kit) given in (|28|) . and the 
corresponding estimated distribution function. Here, the bandwidth appearing in (|2"3"|) has 
been empirically fixed as h — 75, in agreement with the remarks in Application 14.11 We 



choose the following baseline non-monotonic hazard rate function: 



13 



(3.5-lQ- 6 t if < t < 650, 

(29) r(t) = I -4.07143 ■ 10~ 6 1 + 0.00492143 if 650 < t < 1000, 

U-10- 6 t- 0.00715 if t> 1000. 

Fig. [5Jc) shows a realization of the confidence band (|2"6"|) obtained from the above data for 
a = 0.025. For c = 0.00025 the strip r(t) ± c falls within the confidence band, with a 
0.95-level confidence. Condition (|27| is thus fulfilled (see Fig. Eld)). Hence, we are finally 
led to consider model ^ as defensible for the observed data, with r{t) given in (f2"9"| and 
c < 0.00025. 
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